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Abstract 

■ Nambu-Goto action of bosonic string predicts the quark-antiquark potential to be 

Q , V {r) = ^ + ar + ^Q. The coefficient7 = —^^^^4^ is the universal Liischer coefficient of 

\ the Liischer term ^ , which depends upon the space-time dimension 'd'. We take linear 

' term in potential as parent and Liischer term as perturbation for the generation of wave 

^ ^ . function for meson in d space-time dimension. The wave function comes out in terms of 

Airy's infinite polynomial series. With this wave function in higher dimension, we then 
study the Isgur-Wise function for heavy-light mesons and its derivatives. 
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1 Introduction: 

^ I It is well known that in the non-perturbative regime of QCD, i.e, in the low-energy regime, 
00 ' phenomenological models are very useful for the study of hadrons [1]. Since long, potential 
models for mesons, with one heavy and one light quarks, are under theoretical investigation for 
the study of properties of hadrons like their mass, form factors, decay widths etc. 
Q It is worth mentioning here that, for getting the mesonic wave-function in such potential model 
^ approach, the choice of potential is of utmost importance. Several potentials are there for 
• • i modeling quark-antiquark bound states like Martin Potential [2], Cornell Potential [3], Richard- 
• ^ ' son potential [4], Logarithmic potential[5]. Out of these, the Cornell potential, with linear plus 
^ Coulombic form, has been very popular and useful for such phenomenological study in QCD. 
d . The wave function for heavy-light mesons have been calculated earlier with such a potential 
by applying quantum mechanical perturbation technique [6] . This has been deduced both with 
Coulombic term in potential as parent [7-9] and also with linear confinement term as parent [10- 
11]. In the present work, we report wave function for meson developed within some higher 
dimensional string inspired potential model approach. 

In classical string model for hadrons[12-14], proposed by Yoichiro Nambu and Tetsuo Goto, two- 
quark potential can be expressed as [15] V{r) = ^ -\- ar -\- fiQ. Here, the coefficient7 = — ^^^^j^ is 
the universal Liischer coefficient of the Liischer term ^ [16], which depends upon the space-time 
dimension d. a is the string tension which characterizes the strength of the confining force be- 
tween static charges. Its value is 0.178 GeV^ [17,18], /Xq is a regularisation dependent constant. 
This Nambu-Goto potential in higher dimension is analogous to the linear plus Coulombic type 
Cornell potential in 3 -|- 1 dimensional QCD. 

Earlier, we have developed similar formalism considering Liischer term as parent in perturba- 
tion technique [53]. However, physical basis of string potential model supports the fact that the 
linear confinement term as the main contributing part in inter-quark potential and Liischer's 
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term arises in this potential as a first order correction to the hnear confinement term ar [13]. 
In our present approach of developing wave function for meson using perturbation technique, 
we consider the leading linear term as parent and Liischer term as perturbation and use D- 
dimensional Schrodinger equation for solving out the wave function of mesons. Several methods 
for sohition of higher dimensional Schrodinger equation with different potentials are there in 
literaturc[19-29]. This wc take as a generalisation of our previous works in three-dimensional 
QCD [10-11] to higher dimension. 

However, due to the lack of exact solution for D dimensional Schrodinger equation with linear 
plus Coulombic type potential, we opt for employing perturbation method. With linear term 
in potential as parent, the wave function comes out in terms of Airy's function [30], which is 
an infinite polynomial series. Earlier, some approximated wave function for meson has been 
deduced within QCD framework, with linear term in Cornell potential as parent and coulombic 
term as perturbation, considering Airy's polynomial function up to O(r^) [10], which is recently 
further improved upon by the authors considering complete Airy series [52]. 
In this paper, we have reported the wave function for meson in higher space-time dimension con- 
sidering complete Airy's infinite series, following Dalgarno's method of perturbation[6, 31-32]. 
Based on this wave function, we then have studied Isgur-Wise function (IWF) for heavy-light 
mesons [33-34] and its derivatives in higher space-time dimension. IWF is a universal function 
representing all the form factors of heavy-light mesons in the infinite mass limit in semileptonic 
transitions [35] of these mesons. It is a fundamental quantity in QCD, which can be determined 
non-pert urbatively. 

In such a study with wave function containing Airy's infinite series,the infinite upper limit of 
integration in the calculations of normalisation constant and derivatives of IWF, gives rise to 
divergence. In our earlier works [11] we have successfully introduced some reasonable cut-off 
to upper infinite limit of integration in the study of IWF and its derivatives. Here we fix the 
cut-off to infinite upper limit by applying the convergence condition of wave function in per- 
turbation theory. We also study the dimensional dependence of these parameters and compare 
our present results with our earlier ones and also with the corresponding theoretical and ex- 
perimental expectations for D = 3 [36-43,53]. 

With this introduction, detailed formalism is reported in section 2, calculations and results in 
section 3. Final conclusion and remarks are reported in section 4. 



2 Formalism: 

2.1 Potential Model: 

As expressed in equation (1), the quark- antiquark potential for bosonic strings as enunciated 
by Nambu and Goto has the standard form: 

V{r) + + (1) 

with 7 as the Liischer term given by : 

->=-^-^ (2) 

/ 24 ^ ^ 

For the sake of simplification in formalism and comparison with Cornell potential [3] , we take 
Nambu-Goto potential of equation (1) as: 

\ 7 • T '^{d — 2) , . 

V[r) = \- ar + IJ.0 now with 7 = . (3j 
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Hcrc.d is the space-time dimension with d = D + 1, D is the spatial dimension. 
We take err as parent and — ^ as perturbation, with regularisation term = 0. Our unper- 
turbed Hamiltonian is[19]: 

H. = -^^.r (4) 
The Schrodinger equation in D-dimension is [23]: 

[-^V^ + Vo{r)]'^{r, Qn) = ^*(r, Qn) (5) 



with 



and 



^{r,nD) ^ R{r)Y{nn) (7) 

Here, is a generahsation of the centrifugal barrier [44] in D dimension. 

The eigen values of K\,{Q,d) age given by : 

Kl{nD)Y{nD) = i{i + D- 2)Y{nD) (8) 

Here Y{D,d) and R{r) are the spherical harmonics and spherical coordinates; I is the angular 
momentum quantum number and E is the energy Eigen value. 
This gives the equation (5) in term of radial part as : 

Taking h — 1 and for / = state, we have: 



[^ + ^^ + M^-W(r)=0 (10) 

R"{r) + ^^R'{r) + 2fi{E - ar)R{r) = (11) 
r 

We take: 

R{r)^r^U{r) (12) 
Then our equation (11) transforms into : 

^"(0 - ^^'^11?~^^ U{t) + 2ii{E - ar)U(r) = (13) 

or U"ir) - 'I Uir) + 2^xiE - ar)Uir) = with A = (14) 
r"' 2 

Wave function with only linear term in potential : 

We take : 

Q{r) = Q\r - go (15) 
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So, the equation (14) reduces to : 
Case-I: 

For the case ^ — >■ oo , term vanishes, the above equation (16) becomes : 

^M^) = (17) 

Following our previous work [11], we take Qi — {2iiaY/^ and Qq — (^y/^E, so that our 
construction of Q{r) becomes: 



Q{r) = Qir - go ^ {2ixaY'h - {^f'^E (18) 
Then equation (17) transforms into: 

(fU{Q) 



- qU{q) = (19) 



Here, we mention that, are the zeroes of the Airy's function such that Ai[^o] — 0. The 
solution of the equation (19) comes out in terms of Airy's function Ai[Q], as[45]: 



U{q) - Ai\q] = A^A^[(2^(J)^/V - {^f'^E] (20) 



This is the asymptotic solution of equation (17). 
Case-II: 

If we now take p — )■ 0, then the centrifugal term of (g^gp)2 will be dominant so that equation 
(16) can be reduced as : 

A(A+ 1) 

The non-singular solution of this equation is : 



V» - i^,U(d = (21) 



u(e) ~ (e + Qo)-^^^^ (22) 

We construct the approximate analytic solution of equation (16) as the product of two solution 
of extreme cases [46] as: 

U{Q)^{Q+QorAi[Q\ (23) 

The unperturbed ground state meson wave function with linear term in potential, has thus the 
form : 

*°(r,L>) =Arr^(^ir)"^Ai[^ir-^o], m = VA(A + 1) (24) 
Here, N is the normalisation constant for unperturbed wave function. 



Wave function with Liischer term as perturbation: 
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Given this unperturbed wave function in eqn (24), we now proceed to construct the perturbed 
(and hence total) wave function for meson taking Liischer term in potential as perturbation. 
Our D-dimensional Schrodinger equation is: 

[-^(|^ + + - EWir) = (^ + W')nr, D) (25) 

In terms of radial wave function -Ri(r, D), 

[-^(^ + ^^^) + - E)]Ri{r, D)^{^ + W')r^{g^rrA^[g^r - Qo] (26) 



We assume, 



R,{r, D) = r^F{r, D){Q,rrAt[Q,r - g^] (27) 



Employing Dalgarno's method of perturbation the unperturbed wave function comes out as 
(Appendix-A): 

^\r,D) = Nir^{Qir)"'[Ai{r,D)r + A2{r,Dy + A3{r,Dy + ]Ai[gir - go] (28) 

Ai{r) , A2{r), As{r) etc are having explicit form as given by equation (A.24-A.28) of Appendix- 
A. Our the total wave function comes out to be: 

**°*(r,D) = ^°(r)+*'(r) = Nir^[l+Ai{r, D)r+A2{r, Dy+A3{r, 0^+ ]{gir)"'Ai[gir-go] 

(29) 

Here, A^i is the normalisation constant for the total wave function. 



2.2 Isgur-Wise Function and its derivatives: 

For meson containing one heavy quark (c,b,t), the mass of heavy quark is much greater than 
the QCD scale parameter Aqcd', the four-velocity of heavy quark is almost same as the four 
velocity of meson and in the meson rest frame, the heavy quark appears as a static colour 
source[47]. This brings spin flavour symmetry and suppression of form factors[48], resulting in 
a single universal form factor - the Isgur-Wise Function (IWF). This IWF ^{v, v') is a function 
of four velocities v and v' of the heavy particle before and after decay. It is normalised at zero 
recoil. \i y = v.v', then, for zero recoil {y — 1), ^{y) — 1. In terms of this Lorentz boost y, the 
IWF can be expressed as [49] : 

C{y) = 1 - p\y -i) + c(y-ir + (30) 

Here is the slope ( charge radii ) and C is the curvature ( convexity parameter) of IWF, 
which are measured at zero recoil point as : 

2_ say), r^_my)^ 

The calculation of and C provides a measure of the vahdity of HQET along with a vahd test 
for confirmation of our wave function, in infinite quark mass limit. There have been several 
attempts to calculate and C from theory and models[36-43]. The corresponding results are 
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shown in Table-9. The calculation of this IWF is non-perturbative in principle [50]and this 
function depends upon the meson wave function and some kinematic factor, as given below for 
D spatial dimension [54]: 

POO 

/ DCijr^-^|*(r)|^cos(pr)rfr (32) 
Jo 

where Cd^ ,"0 , and cos(pr) = 1 - + ^ + with = 2y[/^(t/ - 1). Taking cos(pr) 

up to O(r^) we get, 

^{y)^DCn r'-'mr)\'dr-[DCDl^' r^+'\^(r)\'dr](y-l)+[^DCDfi' r''+'Mr)\'dr](y-lf 
Jo Jo t) Jq 

(33) 

Comparing (30) and (33) give us : 

POO 

p' = DCD^i^ r^+i|^(r)|2(ir (34) 



^0 

1 f°° 

C = -DCd^i^ r^+='|*(r)Prfr (35) 
6 Jo 

POO 

DCd / r^-^\^{r)\^dr^l (36) 
Jo 

Equation (36) gives the normalization constants A'" and iVi for ^°(r, D) and ^*''*"^(r, D). 

3 Calculations and Results: 

With the wave function constructed, we now proceed to study IWF and its derivatives like slope 
and curvature parameters. Prom expression of total wave function in equation (29), we find that 
it contains two infinite series - one of power series in r with coefficients Ai{r, D), A2{r. D), A^^r, D) 
etc and another of infinite Airy series As the infinite limit of integration in the calcula- 

tions of IWF and its derivatives, involving Airy's infinite series, gives rise to divergence, we opt 
for some reasonable cut-off to this infinite upper limit of integration. In principle, this cut-off 
To should be greater than the size of hadrons ,i.e, ro should be greater than being 
the mass of hadron. The convergence condition of the total wave function obtained through 
perturbation technique implies that ^'(r) < ^(r). From this we obtain: 

I ^i(r, D)r + A2{r, + A^{r, D)r^ + |< 1 (37) 

This condition gives us the limiting values of cut-off ro for different D values. Considering two 
terms in equation (37), we obtain the equation relating cut-off ro and dimension D as: 

A(ro, D)ro + A^iro, D)rl = 1 (38) 

This transforms into : 

K^i{D)tI + K^2{D)tI + Ki3(/?)ro + i^i4 = (39) 

Here, 
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Kn{D) = 2fj,W'f2o (40) 

K,2{D) = -4/^7/10 + 2fj,W'fn - 2//7/20 + '2fj,W'f2i (41) 

Xi3(L') = 2//W7io - 4/x7/n - 2/X7/21 + + 1 (42) 

KuiD) = -2/^7/22 (43) 

f,o{D) = 2/(D) + (D - 1) + 2^>i (44) 

/n(D) = 2^iA: (45) 

f,o{D) = /(i^)[/(i^) - 1] - (D - !)/(!)) + 2fiD) + {D- l)g, (46) 

f2i{D) = {D~l)g^k (47) 

/22(D) = (48) 

m = ^^4^ (49) 

Equation (39) is a cubic equation having only one real root which is given by : 

K12 2'/\-Kf, + 3K^,K,s) , 1 



with 



g = 9KnKi2Kis - 2Kl^ - 27Kl,Ku + ^4(3i^ni^i3 - Kf^f + {9KuKuK,s - 2Kf^ - 27Kl^K^^Y 



Following equation (50) and (40-49), the numerical values of cut-off Tq for D and B mesons for 
different dimensions are shown in Table-1. Also, consider at ion of such cut-off to upper limit 
of integrations will not sacrifice the nature and value of IWF and its derivatives, because, 
Airy's function falls very sharply (almost exponentially) and almost dies out with increas- 
ing r- value [45,51]. In fact, the Airy's function value becomes neghgibly small for r > 4 
{AiryAi[A] — 0.000952). Further, we restrict our calculation up to Airy order and term 
Af,{r,D) in the second infinite series. The values of unperturbed energy E and of perturbed 
energy W for different D value are shown in Table-2. The table reflects that perturbed energy 
W is less than the unperturbed energy i?, as it should be. 

It is to be mentioned here that, Ai{r, D),A2{r, D),A^{r^ D) etc appearing in the second infinite 
series of the wave function are functions of Ci(r) and C2(r). Now, for the simplification of 
calculation, we consider lowest Airy order in computing Ci(r) and C2(r)[Appendix-B]. 
With these more explicit forms of functions Ai{r, D) in meson wave function as in equation (29), 
we now explore IWF and its derivatives for Airy's polynomial order r^*^, considering specific 
cut-off values for specific dimension value D. Results for B and D mesons are shown in Table-3. 
The variation of IWF with lorentz boost y for different D values are shown in Figure- 1. The 
figure is self explanatory, obviously indicating that the zero recoil condition, ^{y) — 1 &X y — 1 
is maintained all-through. 

For comparison of our results for D = the results of different standard models and collabo- 
rations in three dimensional QCD are shown in Table-4. 



4 Conclusion and remarks: 

In this work, we have developed wave function for mesons m d = D + 1 space-time dimension 
considering linear confinement term of Nambu-Goto potential as parent and Liischer term as 
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perturbation. We have then further extended our analysis in finding the derivatives of Isgur- 
Wise function which are here obviously dimension dependent. 

In calculating the slope and curvature of IWF, we have introduced cut-off value to the infi- 
nite upper limit of integrations to overcome divergences; this cut-off is found to be dimension 
dependent. We put forward our observation that our choice of specific cut-off value in these cal- 
culations does not compromise with the accuracy of our results. This will be more established, 
if we consider the asymptotic form of Airy's function (at D = 3) [51]: 



With this asymptotic form (for D = 3)we have also calculated the derivatives of ^(y) consid- 
ering limit of integration from tq to oo. The results for different cut-off values are shown in 
Table-5. The values of and C are exceptionally low here as compared to our calculated 
values thus justifying our consideration of restricting the upper limit of integration to some 
reasonable finite value. 

While studying the compatibility of our results for D — 3 with the standard results of different 
models in 3 dimensional QCD (Table-4), we find that, our results are comparatively lower than 
the expectations. However,while comparing our results with that of our previous work [53] , we 
put forward our comment that in ref[53] the results are exceptionally higher than the expec- 
tations and in our present approach we have overcome this limitation. Although our results 
are comparatively lower than the expectations, we find that and C values gradually go on 
increasing with the increase of dimension D, as is evident from Table-3, which is further con- 
firmed from Figure-2. Further, Figure-1 confirms that the boundary condition of IWF, ^(1) = 1 
is preserved everywhere. 

Lastly, we make the following comments on our present formalism and calculations. 

1. In this work, we have developed meson wave function with an higher dimensional outlook 
and made studies of IWF and its derivatives within such approach. 

2. We opt for the quantum mechanical perturbation technique in deducing the meson wave 
function, due to the constraint in getting exact analytic solution of Schrodinger equation 
involving linear plus Coulombic potential. 

3. As far as numerical accuracy is concerned, the present perturbation approach appears 
short of the numerical solution of Schrodinger equation or calculations of Lattice QCD. 
Still, we believe that the potential model approach gives sufficiently reasonable physical 
insight into the problem. 

4. Lastly, terms Ai{r, D), A2{r, D), A^lr, D) etc. of the second infinite series in the meson 
wave function arc truncated up to A^{r, D) in our calculation. Improvement in the numer- 
ical analysis can be carried out by considering higher terms in both infinite polynomial 
series involved in the wave function. 
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A Appendix 

Determination of equation (28) using Dalgarno's method of perturbation. 

Prom equation (27), we have- 

Ri{r) = gTr'^^'^F{r)Ai[g] (A.l) 

We obtain- 

R'(r) = g^r'-^^[^—^-^^-F(r)Ai[g] + F'(r)Ai[g] + giF(r)Ai'[gm.2) 

2 r 



i-o+2m .1 — D + 2m ,1 — D + 2m _,1 , .r . ^, .. ^1 — D + 2m 1 ^, , , , , 
•^^[ ( l)-Ai[g]F{r) + 2 -F'{r)Az[g] + 

1 — D + 9m 1 

+2 g,F{r)Aif[g\ + F"{r)Ai[g\ + 2g,F'{r)Ai'[g\ + ^?F(r)^z"[^]A.3) 



We take 



We then have- 



Ai'[g\ = Z{r)Ai[g\ so that (A.4) 
Af[g\ = Z''{r)Ai[g\ + Z'{r)Ai[g\ (A.5) 



K{r) = gTr'-^^AM^-^-^-Fir) + F'{r) + g,F{T)Z] (A.6) 

2 r 



i?'/(r) = g^r^^Ai[g][ ( l)^F{r) + 2 -F'{r) 

^ — D + 2m 1 

+2 pF{r)Z + F" + 2g^F\r)Z + glF{r)Z^ + ^?F(r)Z'] (A.7) 

2 r 

Putting equations (A.l), (A.6) and (A.7) in equation(26), we obtain: 



F"{r) + F'{r)Ki{r, D) + F{r)K2{r, D) - 2ii{ar - E)F{r) = 

-2^Ji{- + W') with (A.8) 

K.{r. D) = 2 '-D + ^rn l ^ [^^^ ^ D-l 
2 r r 



1- D + 2mA- D + 2m ,1 1- D + 2ml 

— 2 — ( — 2 1)^+2 — 2 — 

2/^2 . -/X . D -11- D + 2ml , D-l 



+gl{Z^ + Z') + + g,Z (A.IO) 

r 2 r r 
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We take: 



Z(r) = ^^ and (A.ll) 
Z\r) + Z'{r)^^ (A.12) 

Our Ki{r, D) and K2(r, D) transform to : 

1 — D + Im 1 1 

K,{r, D) = [2 + {D-l) + 2Q,C,{r)]- = M,{r, D)- (A.13) 

+ (2 ^' + (I? - 1))qiC,{t) + ^?C2(r)]l = M2(r, D)-^ (A.14) 

where- 

Mi(r, n) = 2 ^^-il^ + (£1 - 1) + 2e,Ci(r) (A.16) 

1 — D -I- 9777 

+(2 + (L> - l))^iCi(r) + QlC^ir) (A.16) 

It is to be noted that at D = 3, 

i^i(r)=2^iCi(r)- (A.17) 
r 

i^2(r) = ^?C2(r)i (A.18) 

And equation (A. 8) becomes: 
F"{r) + 2^iCi(r)iF'(r) + ^?C2(r)^F(r) - 2//((7r - E)F{r) = -2//(^ + W) (A.19) 
which is similar to the equation obtained with 3-dimensional QCD potential. We take: 

F(r, D) = Y^ Ai{r, Oy (A.20) 

so that - 

F'(r, D) = lY, Mr, Dy~' (A.21) 
I 

F"{r, D) = l{l - 1) 5^ Mr, Oy-^ (A.22) 
I 

Applying equations (A.20-A.22) in equation (A.19), we obtain: 

Ai{r, D)[l{l-1)+Mi{r, D)l+M2{r, D)y-^+2iiEEiAi{r, D)r^-2//a ^ Ai{r, Dy+^ = -^-2//W^' 

(A.23) 
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Equating power of r ^ on both sides of equation (A. 23), we get A0M2 — 0, which imply that 
^0 = 0. 

Further, equating powers of ^r^ ^r^r^ and of equation (A. 23), we obtain : 

^) = (6 + 3M, + M.)(M, + M.) ^^-26) 
V^^XM + M2) - Va7(2 + 2Mi + M2) 

^ (12 + 4Mi + M2)(2 + 2Mi + M2)(Mi + M2) ^ ' ' 
^ 8fi^E^j{2 + 2Mi + M2) + 4/i^7iyX6 + 3M, + M2)(Mi + M2) 

(20 + 5Mi + M2)(6 + 3Mi + M2)(2 + 2Mi + M2)(Mi + M2) ^ ' ^ 

With these expressions of Ai{r, D),A2{r, D),As{r, D) etc we can now construct F(r, D) as: 

F(r, D) = Ai{r, D)r + A2(r, D)r2 + ^(r, D)r2 + A4(r, D)r* + ^(r, D)r^ + (A.29) 

And our radial wave function i?i(r, D) comes out to be: 

Ri{r,D) = r^{gir)"'F{r,D)Ai[g] (A.30) 

The perturbed wave function ^'(r) and total wave function ^*°*(r) are constructed as below: 

*'(r> D) = Nir^ (^ir)'"[A(r, D)r + A2{r, Dy + A^^r, + ]Ai[Qir - ^oIA.31) 

As, ¥°\r,D) = ^>^{r,D) + m'{r,D) 
= AAir^(^ir)™[l + Ai{r, D)r + ^2(r, D)r^ + ^(r, + ]^^[^lr - ^oIA.32) 

B Appendix 

Calculation of Ci{r) and C2(r) ; 
The Airy's infinite series as a function of ^ = gir — go can be expressed as [51] : 

- = + 6 + 180 + 12960 + - 

^ , iQir-Qo)'^ , (gir - gp)^ , (gir-go)lO 
^^[^^^^ - + 12 + 504 + 45360 + ^^"'^ 

with oi = 32y3r(2/3) = 0.3550281 and h = ^TTs^^Yy^y = 0.2588194. 
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Table 1: 


Values o: 


To for different D 


D meson 


B meson 


D 


ro 


D 




3 


3.6354 


3 


3.4719 


4 


3.9482 


4 


3.7678 


5 


4.1361 


5 


3.9437 


9 


4.6362 


9 


4.4135 


15 


5.1158 


15 


4.8657 


20 


5.4104 


20 


5.1436 


25 


5.6515 


25 


5.3712 



To find Ci(r) and C2(r) , we take truncated Airy series up to lowest order, so that, we have : 

Z(r) = ^ = -^^^(^^^ - (B.2) 
r Ai{g,r -go) 



Also, 



We thus obtain- 



-biQi 



oi - bi{gir - Qo) 
biQi 

bi{gir - Qo) - ai 



(B.3) 
(B.4) 



= ill-^l±^V (B.5) 
r biQi r 

= 1[1-V^ = -(1 + -) (B.6) 

rj^ ry ry ly* 

k 

Therefore, Ci(r) = 1 + - (B.7) 

with k = ^l±h^ (B.8) 
biQi 



1 2A; 

an(i Z'(r) = — so that, (B.IO) 

^-Z\r)^Z'[r)^^^^ (B.ll) 



C2{r) = ^ (B.12) 
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1.6- 
1.4- 
1.2- 
1.0- 
J- 0.8- 
0.6- 
0.4- 
0.2- 



- D meson 
B meson 




0.7 

0.6 

0.5 

g) 0.4 - 
3 

13 

£ 0.3 

o 

0.2 
0.1 
0.0 



- D meson 
B meson 



(a) vs D for D meson 



(b) p vs D for B meson 



Figure 2: Variation of and C with D for D and B mesons. 
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Tabl 


e 2: Values of E and 


W in 


GeV for different D 


D meson 


B meson 


D 


E 


W 


D 


E 


W 


3 


0.345938 


0.137536 


3 


0.352764 


0.139463 


4 


0.637435 


0.141397 


4 


0.708288 


0.143216 


5 


0.832819 


0.148987 


5 


1.0102 


0.149692 


9 


1.54202 


0.163230 


9 


1.66759 


0.16872 


15 


1.73322 


0.195488 


15 


1.74195 


0.20805 


20 


1.7477 


0.240184 


20 


1.74838 


0.24205 


25 


1.75193 


0.306038 


25 


1.76108 


0.306265 



Table 3: and C for different D values. 





D meson 


B meson 


D 


p2 C 


p2 C 


3 


0.2158 0.0174 


0.2608 0.0254 


4 


0.3874 0.0457 


0.4039 0.0483 


5 


0.1366 0.0613 


0.1813 0.0679 


9 


0.8355 0.1746 


0.9497 0.1916 


15 


0.9622 0.2710 


1.0281 0.3274 


20 


1.2178 0.3997 


1.3352 0.4887 


25 


1.3732 0.5041 


1.4637 0.6623 



Table 4: Results of slope and curvature of ^{y) in different models and collaborations. 



Model / collaboration 


Value of slope 


Value of curvature 


Ref [10] 


0.7936 


0.0008 


Le Youanc et al [36] 


> 0.75 


> 0.47 


Skryme Model [37] 


1.3 


0.85 


Neubert [38] 


0.82 




UK QCD Collab. [39] 


0.83 




CLEO [40] 


1.67 




BELLE [41] 


1.35 




HFAG [42] 


1.17 ±0.05 




Huang [43] 


1.35 ±0.12 





Table 5: Values of and C with asymptotic form of Airy's function. 



To value 


(asymptotic ) 


C (asymptotic) 


5 


4.6 X 10-9 


1.6 X 10-*^ 


6 


5.027 X 10-" 


2.464 X 10-" 


7 


3.56 X 10-^=^ 


2.345 X 10-^3 


8 


1.695 X 10"^^ 


7.028 X 10"^^ 


9 


5.248 X 10-^« 


6.597 X 10"^^ 


10 


2.92 X 10-^9 


2.78 X 10-^9 
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